Simultaneous decompositions of two states 



O 

o 

^ ■ Armin Uhlmann, Leipzig [| 



> 
oo 

o 

o 



X 



Dedicated to Roman S. Ingarden 



Abstract 



Simultaneous decompositions of a pair of states into pure ones 
. are examined. There are privileged decompositions which are distin- 

I guished from all the other ones. 

, Presently we witness that quantum information theory is becoming an inter- 

qh| disciplinary, quickly growing field of research. In its history Roman S. In- 

garden is playing a significant role, both by his own research and by posing 
stimulating questions and problems [|I|. It is about 40 years ago that I met 
Roman the first time, and he was already thinking about the role of infor- 
mation in quantum physics and, in particular, whether one can found the 
concept of probability onto that of information [0]. 
I feel honored by the possibility to dedicate to him the following paper. 



1 Decomposing one density operator 

A density operator, representing a state, is a positive operator with trace 
one. However, it is convenient for the following considerations not to insist 
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in normalization. 

We shall assume, mainly for technical simplicity, a finite dimensional Hilbert 
space, ?i, the dimension of which is denoted by dim?i = d. Thus, math- 
ematically, we are just dealing with positive operators (and with the Null 
operator) of a finite dimensional Hilbert space. 

Let r be a positive operator on our Hilbert space. Its decreasingly ordered 
eigenvalues are denoted by Ai, A2, . . ., i. e. 

spec(r) = {Ai > A2 > . . .} 

By a decomposition of t I denote every set of vectors \xj) such that 

^ = Elx.)(x.l (1) 

As I showed in 

m m 

EA.>E(x.lx.) (2) 

is valid for all 1 < m < dimTY. Moreover, equality is reached if and only if 
\Xj) is an eigenvector for Xj of r for all j = 1, . . . , m. 

The motivation for asking questions of that kind has been the problem 
whether the von Neumann entropy of a density operator is already fixed 
by its position as a point in the convex set of all density operators. The 
result just quoted gives, if written with normalized vectors, an affirmative 
answer. Indeed, my aim was to define on every (compact) convex set a func- 
tion which just gives the von Neumann entropy if applied to state spaces of 
a quantum system. Up to day I do not know whether the construction is of 
any use for other convex sets than quantum state spaces. 

In M. A. Nielsen proved the reversed statement: If are positive numbers 
which are majorized by spec(T), then there exists a decomposition (|l|) such 
that pj = ixjlXj)- 

The results mentioned above will be slightly extended to the case that there 
are two decompositions of one and the same r. Thus let 

^ = ElXj)(Xjl (3) 

be a further decomposition of r. Adding (|1|) and (|^) we get a decomposition 
of 2r and 

m m 

E((x.ix.) + (x;-ix;))>2Ei(x.ix;-)i 
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Equality takes place iff \xj) differs from by a phase factor only. Because 
the eigenvalues of 2r are just 2Xj we get 

Proposition 1: Let and ^ he decompositions of t and Ai > A2 > . . . 

the decreasingly ordered eigenvalues of r, and 1 < m < d. Then 

m m 

E^.>Ei(x.ix;)i (4) 

j=i i=i 
Equality holds if and only if for 1 < j < m 

^\Xj) = ^j\Xj), \x'j)=^j\Xj) (5) 
with unimodular numbers ej. 

2 Decomposing two density operators 

Let us now consider a pair, p and u, of positive operators. 
Definition: F^{p,u) denotes the sum of the m largest eigenvalues of 

(v^^v^)^/' (6) 

The definition works well for 1 < m < d. It is sometimes convenient to 
extend it by = ii m > d and to set = for m = 0. 
Remark that is the square root of the transition probability The 
square root of the transition probability is called fidelity and is denoted by 
F{p, u) in the present paper. Notice, however, that Jozsa, who showed its 
use in quantum information theory 0, identified the general transition prob- 
ability with his fidelity concept (and not with its square root). 
A further remark is the following: In [0] I considered another quantity: The 
fc-fidelity, Fk, which is the sum of all but the first k eigenvalues of (H). These 
partial fidehties are jointly concave (and super-additive) in its arguments for 
k = 0, 1, ... Obviously, 

F+{p,uj) = F{p,uj)-FUp,io) 

In contrast to the partial fidelities, the quantity (|[) seems to be neither 
concave nor convex if m is smaller than dimTi. 
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Finally, let us rewrite (^) of proposition 1 as 

m 

nr,r)>Y.\{xM)l (7) 
i=i 

Remember that equality in can be reached by eigenvector decompositions 
of r with decreasingly ordered eigenvalues. 

Theorem 1: Let be 1 < m < d. It is 

m 

F+(p,a;) = max^|(^,|^,)| (8) 
i=i 

where the maximum is to perform over all possible decompositions 

// the length of a decomposition is less than dim Ti, or if the length of the two 
compositions are different, one adds some zero vectors to get decomposi- 
tions of equal and large enough length. 

The proof of the theorem starts by stating the invariance of the eigenvalues 
of with respect of a transformation 

{p,u} {p,a;}^ := {XpX*,{X-'yujX-'} (10) 

for any invertible operator X, see |0. ((In the present paper the Hermitian 
adjoint of an operator A is denoted by A* and not by A^.)) Hence the sum 
of the m largest eigenvalues of (^) cannot be changed by such a transforma- 
tion. On the other hand, if we simultaneously transform decompositions (^ 
according to 

\^^) ^ (x-'y\^) (11) 

then the right hand side of (P) remains unchanged. Therefore, if the assertion 
of the theorem is true for a pair of density operators {p, u}, it is true for every 
pair {p, u}^ . 

Let now p and uj be invertible (i. e. faithful). If we then can choose X such 
that 

{p,a;}^:={r,r} (12) 
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with a certain r yet to be determined, we are done: For the pair {r, r} the 
theorem is equivalent to proposition 1. But 
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:= T 



(13) 




(14) 



Hence, the theorem is true for invertible p and u. 

Indeed, the proof covers the case of any pair p, u, with equal supports: To 
see it we only have to replace H by the supporting Hilbert subspace because 
neither to nor to the decompositions there is a non-zero contribution 
from the null spaces (i.e. the kernels) of p and u. 

We now prove that the right hand side of (P) never exceeds F+. Denote by 
Po, Qo the projection operators onto the null spaces of p and u. We choose 
decompositions of Pq and Qo with vectors lipi) and respectively. We 
complement arbitrarily chosen decompositions (|^) to those of p' = p + CiPq 
and uj' = u + C2Qq with cj > 0. For p' this is done by 



and similarly we proceed with u'. Because p' and u' are invertible, we already 
can apply theorem 1 to them. Because F^{p',u') is approaching F^{p,u) if 
Cj — > we are done. 

What remains to show is the following: There are decompositions @ such 
that J2 I ('04: \ is equal to F+, whatsoever the support properties of p and 
uj may be. To get this we first assert: 

Let Q be the projection operator onto the supporting space of u>. For all 
decompositions (H) we get 



because every vector of a decomposition of u must be an eigenvector of Q. 
That is, every one of the sums in question for p,u gives one for QpQ,u} 
yielding the same value. On the other hand, if we start with decompositions 
of QpQ, u, we can add terms orthogonal to tu into the decomposition of p to 
get a decompositions of p, uj without changing the value of the sum. Below 
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we shall show the equality of F^{p,u) with F^{QpQ,ll!), and, all together, 
we obtain: If and only if theorem 1 is true for the pair QpQ, uj it is true 
for the pair p, oj. Now we can proceed as following: If the supports of QpQ 
and OJ are equal, we are done. If not, we consider the projection operator 
Pi onto the support of QpQ, yielding the same statement for the pairs p, a;, 
QpQ,u, and QpQ, PituPi. Either the last pair is of equal support, and we 
are done, or continue the same game with the projection operator Qi onto 
the support of Pi a; Pi. This procedure must terminate after a finite number 
of steps yielding a pair with equal supports. The obvious reason: In every 
necessary step, the rank of one member of the pair under consideration is 
diminished, and we are in finite dimensions. 

The proof of theorem 1 is done after showing the equality of F^{p,uj) with 
F^{QpQ, u) if Q is the support projection of uj. This assertion is a particular 
case with X = Q of the equation 

P+(p,X*.;X) = P+(XpX*,a;) (15) 

For invertible X the transformation (|lO|) does not change the eigenvalues of 
(|^). By the replacement u — X*ujX we thus get (^) for invertible X. But 
F^ is continuous in its arguments, and (p!5[) is valid for all X. 

Let us underline the main point in constructing decompositions (^ satisfying 

m 

^^(P,^) = E(^.l¥'.), m = l,2,... (16) 

3=1 

We have to solve ([T3|) so that X and r are at our disposal. From the spectral 
decomposition of r, 

7" = m Xi) {Xj I , {Xj I Xk) = XjSjk (17) 
we get an optimal decomposition satisfying (|r^) by 

\^P,)=X-'\x,), \ip,) = X*\x,) (18) 
Such a choice fulfills the bi-orthogonal relations 

i'lpkl'^j) = {ipj\(pk) = XjSjk (19) 

Acknowledgement. I like to thank P. M. Alberti, B. Crell, Ch. Fuchs, and 
M. Nielsen for valuable discussions and correspondence. 
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